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2. 


1.  Introduction.  Let  X.  = (X.  ...... X .)  , iil  be  a 

-l  li  pi 

sequence  of  independent  random  vectors  (i.  rvs)  with  continuous 

cumulative  distribution  functions  (cdfs) 

(1.1)  F (x)  = PC Xi  < x]  = F ( x - a - pCj ) , i i 1 , x f RP  , p > 1 

t 

where  a = (a. , . . . ,a  ) , 0 = ( (0  ) ) . , , q?l  are 

k=l, . . . ,q 
/ 

unknown  parameters  and  c.  = (c,  ......  c .)  , i?l  are  known 

^ -i  li  qi 

vectors  of  regression  constants.  We  partition 

(1.2)  9 = (0^,  £2  ) , qi  + q2  = q * qi  ? 0 ' i = 1*2  • 

PxqL  pxq2 

The  problem  is  to  test  the  null  (composite)  hypothesis 

(1.3)  Hq  : 02  = 0 against  : p2  ? 0 . 

We  may  mention  that  by  the  classical  canonical  reduction 
fviz.  Anderson  (1958),  Chapter  8)],  a general  linear  hypothesis 
on  0 can  always  be  reduced  to  a form  similar  to  (1.3).  For  a 
particular  case  of  q^  = q i.e.  HQ  : £ = 0 , the  problem  reduces 
to  that  of  testing  a simple  null  hypothesis,  the  rank  order  tests 
for  which  have  already  been  studied  by  Puri  and  Sen  (1969) . 
However,  the  technique  developed  in  that  paper  is  not  applicable 
when  q2 < q . This  difficulty  is  circumvented  here  by  using 


3. 


aligned  rank  order  tests  fas  in  Sen  (1969),  and,  Puri  and 
Sen  (1973)  both  dealing  with  the  univariate  models]  where 
the  alignment  is  based  on  estimates  of  £ developed  in  Sen 
and  Puri  (1969)  and  Jurettkova'  (1971). 

The  proposed  rank  order  tests  for  are  considered 

in  section  3 following  the  preliminary  notions  and  basic 
assumptions  in  section  2.  Section  4 deals  with  asymptotic 
comparison  of  parametric  and  rank  order  tests,  and  the 
asymptotic  optimality  of  the  proposed  tests.  The  last 
section  deals  with  a special  case  of  (1.3),  namely,  testing 
the  hypothesis  of  parallelism  of  several  multiple  regression 
surfaces  which  turn  out  to  be  the  multivariate  multi- 
parameter analogue  of  Sen  (1969) . 


n 


2.  Notations  and  assumptions.  Let  R..  =7  u (X . . -X.  ) , 


3i 


a=l 


3i  3o 


(where  u(t)  =1  or  0 according  as  t is  ^ or  < 0)  be 

the  rank  of  X . . among  X ,X  . ; i=l,...,n;j  = l p 

ji  3l  3n 

Since  F is  continuous,  ties  among  the  observations  may  be 


neglected  in  probability.  For  each  j (=l,...,p)  , consider 
a set  of  scores  a^Nl),...,  a^\n)  , generated  by  a function 

ePj(u)  , 0 < u < 1 , in  either  of.  the  following  ways. 


4. 


I 


(2.1)  a ^ (i)  = cp  . (i/(n+l) ) or  a ^ (i)  = Ecp  . (U  .)  , 1 < i < n j 1 -c  j < p 

n j n ] ni 

where  cp  (u)  is  assumed  to  be  square  integrable  inside  (0,1)  , 

and  u < ...  <u  is  an  order  statistic  of  a sample  of  size 
nl  nn 

n from  the  rectangular  distribution  over  (0,1).  Our  proposed 
procedure  is  based  on  the  following  type  of  rank  order 
statistics . 


n 


(2.2)  S = ((S  )) 

n n,  jk 


3 " *I,(C 


(j) 


n,  jk 


i=l 


(R.  .) 


where 


c 


kn 


-1 

n 


n 


k = 1. 


.q  * j 1 , • • • . p • 


Following  Ha'jek  (1968)  and  Hoeffding  (1973),  we  assume 
that  for  every  j(*l,...,p)  , 

(2.3)  q>j(u)  = cp^  (u)  - cp,!2)  (u) 


. (s) 

where  cp^  (u) 

decreasing  in 


, s = l,2  is  absolutely  continuous  and  non- 
up  (0,1)  and 


(2.4)  f |cp(s)  (u)  ] {u(l-u)  }_i5du  <®  • s = 1 , 2 ; j = 1,  . . . ,p 
• 0 D 


Denote 

- r1 

(2.5)  cp  . = | cp.(u)du  , j = l p , 

3 v 0 3 


and 


5. 


00  00 

(2.6)  X..,(F)  -J  i^j(FtJ],X,,tpj'(FI:3'](y,,aP[jj-]t='.y)  -Vfy 

where  Fj._.^(x)  and  F[jj/]^X,y^  are  the  cdfs  of 

jth  and  (j.j')  the  components  respectively.  Assume 


(2.7)  _A  (F)  = ( (X  j j , (F)  ) ) is  positive  definite  and  finite. 


-»  C as  n 


Regarding  the  c^  , we  assume  that 

(2.8)  n V (c.  -c  ) (c.  -c  ) = n_1C 

.L.  ~i  ~n  ~i  ~n  ~n 


where  for  every  n > n^  , 


(2.9)  Cn  = ( (C^  ^/))  is  positive  definite  and  finite 


(2.10) 


„ _ (D  (2) 

~ ~i  ' l**.»*n 


where  for  each  k(=l,...,q)  and  s(=l,2)  . c,.  is  non- 

ki 

decreasing  in  i . (Note  that  the  assumption  (2.7)  is 
a slightly  simplified  version  of  a parallel  assumption  made  by 
Jure^kova'  (1971).  For  q=l  , this  assumption  is  not  necessary) 
Finally,  we  assume  that  for  every  e >0  , there  exists  an 
integer  nQ  = ng(e)  such  that  for  n > nQ  , 

<2-U)  “"^n.lOc’'  f,"ax  lcki-aiJ21  ' 11  = 1 1 • 

l<i<n 


6. 


Regarding  the  cdf  F , we  assume  that  for  each 
j(=l,...,p)  , the  marginal  cdf  ^as  an  absolutely  continuous 

density  function  fj._.^(x)  with  a finite  Fisher  information 

00 

(2.12)  1 j = 1 j ( f j)  =J  {d/dx)log  f (x)}  dF^  (x)  , j=l 

“00 

To  explain  the  alignment  procedure,  we  need  the  following 
notations. 

Let  B = ( (b ..) ) be  a pxq  matrix  with  real  elements 


/ / / 
i = (b  , . . . ,b  ) 

~1  ~p 

(b.)  ) , 1 * i 
Da  ~l 

among  X . (b  . ) 

Da  ~J 


Now  replace  the  R..  in  (2.2)  by  R . . (b . ) for  Ui<n  , 

Ji  Ji  ~1 

1 < j < p and  denote  the  corresponding  matrix  of  rank  order 
statistics  by 


I •••  l 


P • 


and  let 

(2.13) 

X.  (B)  . 
— i — 

X.  - BC.  ; i = l,...,n 
-1  —1 

(2.14) 

R.  . (B) 
11  ~ 

n 

= R.  . (b.)  = V u (X  (b  , ) 

ail  1,1  ^ 

so  that 

VS' 

is  the  rank  of  X . . (b 

li  ~ 

a = 1, . . 

. , n , 1 < i < n . 

(2.15)  S (B)  = («  (b.)))  , j = l,...,p  ; k = 1, q . 

— n ~ n,  jk  —j 

Note  that  by  varying  B on  R^xC*  , we  obtain  a multi- 
parameter multidimensional  stochastic  process  which  is  used 
in  the  next  section  to  introduce  the  proposed  aligned  rank 


order  statistics. 


7. 


3.  The  Proposed  Aligned  Rank  Order  Tests.  As  in  (1.2),  we 
partition  B as 


(3.1)  B = (j^**^)  , Bi  is  pxq.  ; i=l,2  ; qi  + = q 

(3‘2)  Si  = (-^i  (1)  ' Si(2))  ' -Si  (s)  is  a qs-vector,  s = l,2 


Then,  under  HQ  in  (1.3),  we  have 


(3.3)  F^x)  = F (*  ~ ®ifi  (1 ) ) ' 1 < i < n 


First,  we  proceed  to  estimate  jj  for  the  model  (3.3). 

For  this,  consider  the  p x matrix 

(3.4)  S ...  (B  ) = ((S  (b!15)))  . , 

~n(l)  -1  n,  jk  —j  j = 1,  . . . ,p  ; k = 1, q 

where 

(1)  > (2)  ' 

(3.5)  b'  = (b_.  , bj  ) is  a partition  of  b_.  by  (3.1) 

Now  under  (3.3),  has  expectation  o , and  dispersion 

matrix 


(3.6)  A (F)s>C  where  C = 

— ~n(ll)  ~n 


~n(ll)  ' ~n ( 12 ) 1 
~n(21)  ' ?n(22 )l 


(and  ® stands  for  the  Kronecker  product)  and  from  the  results 
of  Puri  and  Sen  (1969),  it  follows  that  for  large  n , under 
the  assumptions  of  section  2, 


8. 


(3.7)  t(n''5Sn(1)  -11)  - npxqi(2-  *<F>  ®£(11)> 

where  £(11)  is  the  qi  x qi  minor  of  c defined  in  (2.8). 
Consequently,  by  the  same  alignment  procedure  as  in  Sen  and 
Puri  (1969)  and  JureSkova'  (1971),  we  define 

P qi 

(3-8)  £n  = {~1  : I I 'Sn,  jk(~j1)}  I = minimum}  . 

j=l  k=l 

Our  proposed  estimator  of  0 (under  (3.3))  is  then 

(3.9)  B = center  of  gravity  of  D 
* n n 

By  arguments  parallel  to  those  of  Jure£kova'  (1971),  it 


follows 

that 

(3.10) 

3UP  ||0  - B ||  g 

B ( D 1,n 

— 1 ~n 

0 t 9S  n -4  co 

(3.11) 

Zin  [^l,n  ‘ -&1])  ■* 

npxq(0'  X(F)  ®£ 

where 

(3.12) 

T(F)  = ( (x j j ' (F) ) ) 

= ( (X  . . , (f)/a  .a 
JJ  j : 

and 

00 

(x)  (x) 

)aFtji  • 

(3.13) 

AJ  =i,d/ax)®3(Pt3] 

9. 


t 


B.  i9  a translation-invariant,  robust,  consistent  and 
asymptotically  normally  distributed  estimator  of  when 

(3.3)  holds.  Our  proposed  tests  are  based  on  the  aligned  rank 
order  statistics 


(3.14)  S = ( (S  ) ) . , 

n(2)  n,jk  ]sl,...iP  ’ k = ^ + 1, ,q 

where 

* $ - (i ) * 

(3'15)  Sn.j*  =i|1(clti-ck.n)an  <V'  ’l*1'1'’' 


(3.16)  R.  . = R.,  (0 


ji  = (£x  n,  0)  , 1 < i <n  , lsjsp 


To  introduce  the  proposed  test  statistics,  we  first  define 


(3.17) 


M = ( (m  , ) ) where 

~n  33  #n 


(3.18)  m , = (n-ll^jf^’lR  )a*]/)(R  ( )-a(i)a‘(j,)} 

jD  # ^ ^ j ^ J i n n J 


where 

(3.19) 


j.y  = i. 


;'j)  . U),p 

" lii 

A 

Also,  replacing  R by  R.,  , l«i<w  , 1 s j - p in 


ji 


ji 


(3.18),  we  denote  the  corresponding  matrix  by 


(3.20) 


j / _) ) • 

~n  3 3 / n 


Let  then 


10. 


(3-21)  ?n  ~n(22)  _ ~n(21)~n(ll)~n(12) 


(3.22)  G = M ®C 
— n — n —n 


pq2  x pq2 


(3-23)  -«„*  "Vjk  P , 

P<32  X bq2 

Our  proposed  test  statistic  is 


(3.24) 


£ = TrfH  G *]  ; 

n ~n— n 


In  the  remainder  of  the  section,  we  show  that  under  H 

0 

in  (1.3)  and  the  assumptions  of  section  2,  has  asymptotically 

a chi  square  distribution  with  pq2  degrees  of  freedom.  This 
provides  an  ADF  (asymptotically  distribution  free)  test  for 


Lemma  3.1.  Under  the  assumptions  of  section  2.  when  It  holds. 


(3.25) 


1-1  p -1 *-l 
nG  5 A (F)  ® C , as  n -* 

~n  — — 


(3.26) 


* -1 

C=C  -C  c c. 

~ (22)  ~ (21)~  ( 11)~  (12) 


★ p * 

Proof.  By  virtue  of  (2.8),  5 C , as  n -*  <*>  . Thus  to  prove 


(3.25),  it  suffices  to  s how  that 


11. 


(3.27) 


Mn  S a (F1)  » as  n -* 


Also  since  m . . , = m . . = (n-1)  V Ta  (i)  l2} 

D 3 .n  D3,n  n n J 

-♦  X . . (F)  = X..  by  (2.1)  and  some  routine  computations,  wo 
33  33 

need  only  to  show  that  for  every  j ^ j ' , 


(3.28)  m..,  5 X • • / (F)  when  holds  . 

33  ,n  0 

By  assumption  (2.3),  (see  also  Ha'jek  (1968),  section  5) 
for  every  e > 0 , there  exists  a decomposition 

(3.29)  cp j (u)  = cp ^ 1 ^ (u)  + cp_!2^  (u)  - cp^3)  (u)  , 0 < u < l , 

3 1,...,, 

. (1)  . , . . (2)  . (3) 

where  cp_.  is  a polynomial,  cp_.  and  cp  . are  non- 
decreasing, and 


3 1 


(3.30) 


V f (u)  ]2du  < e X • • * 1 < j ^ p 

&2  J0  3 33 


Using  (3.29)  we  decompose  n into  9 terms.  Using  the 

Cauchy-Schwarz  inequality  for  the  eight  terms  for  which  at 
least  one  factor  is  non  polynomial  along  with  (3.30),  it  follows 
that  to  prove  (3.28),  it  suffices  to  take  cp^  = cp  ^ , 1 < j «■  p . 
Since  the  cp ^ ^ ^ are  absolutely  continuous  and  are  polynomials, 

A 

for  them,  the  corresponding  m..,  can  be  written  as 

33  »n 


12. 


?<r(l)‘K  (l)''  ~ * 

! cp!  (H  . (x)  )cp  ' (H  , (y)  )dH  ..  (x,y))+o(l)  whe 
•>  »'  J nJ  D nj  njj. 


re 


— CO  —CD 


Hn_.  is  the  sample  cdf  for  the  aligned  observations  on  the  jth 

* * 

variate,  1 < j < p , and  H is  the  bivariate  sample  cdf  for 

njj 

these  aligned  observations.  By  (2.6),  (2.14)  and  (3.11),  on 

ft 

denoting  by  , H^j,  the  corresponding  sample  cdfs  for 

A it  ^ + 

B , it  follows  that  suplH  . - H . . .1  -♦  0 , as  n -♦  ® . 

~1*n  x ' "J!  n,  jj'1 

(i)  ,y 

Also  note  that  the  cp_.  's  are  bounded,  continuous  functions. 

So  first  replacing  H by  H , H by  H , and  then  using 

n n n n 

theorem  4.1  of  Puri  and  Sen  (1969),  the  desired  result  follows. 
In  fact,  it  can  be  shown  that  (3.27)  holds  almost  surely. 


Lemma  3.2.  Under  the  assumptions  of  section  2,  when  H holds. 


(3.31) 


n^ts 


Sn(2,  ' ®n(l)(»l'  S>  +i<B1,n-£l>Sn(12)1  S ° 


as  n ao  , where 


(3.32) 


A = Diag  (A  , . . . , A ) . 

- 1 p 


The  proof  follows  as  a direct  multivariate  extension  of 
Theorem  3.1  of  Jurefckova'  (1971),  and  hence,  the  details  are 


omitted. 


13. 


A L 

By  noting  that  (Bx  n>0)  = Op(n7).(see  Jurefrkova'  (1971)), 

the  following  lemma  also  follow  directly  as  a multivariate 
extension  of  Theorem  3.1  of  Jurefckova'  (1971). 


Lemma  3.3.  Under  the  assumptions  of  section  2,  when  HQ  holds. 


(3‘33)  "^nlll  - Si>2„(ll,]  s 0 • 


as  n -» os  . 

Using  Lemmas  3.2  and  3.3,  we  arrive  at  the  following  result. 


Lemma  3.4.  Under  HQ  in  (1.3)  and  the  assumptions  of  section  2. 


(3-34)  “J,r5n(2)  " -nW  (£l'2>  +Sn(1)(61'2>5MU)Sn(12)lD-  0 


as  n -♦  oo  . 


* 

Consider  now  : £ = 0 . Then  under  HQ  : - 0 , the 

statistics  rs  ,_.(B,#0),  S /,v(B,#°)]  have  the  same  -joint 
— n \ z ) — 1 — —mi;  —1  — 

★ 

distribution  as  that  of  S under  H , and  since  the  later  is 

-n  0 

asymptotically  multi-normal  with  mean  vector  0 and  dispersion 
matrix 


(3.35)  A(P)  ®Cn  , 

it  follows  that  under  in  (1.3), 


14. 


(3.36)  I (n-‘iCSn(2)  (6l,0)  -Sn(1)(g1.0)C^11|Cn(12)]) 

- "pxq2  -•  i(F'*S(22)  -£,2»£Tu)S(1«4 5'  ' 

Hence  using  (3.34)  and  (3.36),  under  in  (1.3),  we  find 

that 


(3.37)  ) -*  h (0,  A(F)®C*) 

n(2)  Pxq2  ~ - ~ 

From  Lemma  3.1,  (3.37)  and  the  asymptotic  distribution  of 

quadratic  forms  associated  with  asymptotically  multinormal 
vectors,  it  follows  that  (under  HQ  in  (1.3)  and  the  conditions 
of  section  2) , 


(3.38) 


Z (X„)  -♦  , as  n -*  <*> 

N pq 


Thus  the  proposed  ADF  test  is  as  follows: 


0 

*0 


if 

Z„ 

v 2 
* X 

N 

P q2»a 

if 

£„ 

<x2 

N 

pq2.a 

where  T is  the  upper  100a%  point  of  the  chi  square 
t,  a 

distribution  with  t degrees  of  freedom. 


4.  Asymptotic  comparison  with  parametric  test.  Consider  now 

a sequence  {Kn}  of  Pitman-type  alternative  hypotheses,  viz. 


15. 


# 

( n) 

(4.1)  K : g_  = B ' = n y0  , y_  is  fixed  and  non-null. 

n — z z ~ z ~z 

Our  aim  is  to  make  the  asymptotic  power  comparisons 
between  the  proposed  rank  order  test  and  the  normal  theory 
likelihood  ratio  test  when  the  underlying  cdf  is  not  necessarily 
multinormal.  Proceeding  as  in  Sen  and  Puri  (1970) . (Where  the 
distribution  theory  of  the  normal  theory  likelihood  ratio 
test  for  the  general  linear  hypotheses  is  considered) , it 
follows  that  if  F possesses  a finite  second  order  moments,  then 
(i)  under  , the  normal  theory  likelihood  ratio  statistic 

i ractually-2  log  (likelihood  ratio  statistic)],  denoted  by 

has  asymptotically  a chi  square  distribution  with  pq^ 
degrees  of  freedom,  and  (ii)  under  {K^}  , it  has  asymptotically 
a non  central  chi  square  distribution  with  pq2  degrees  of 
freedom  and  non-centrality  parameter 

(4.2)  flL  = Tr[f  • (V  (F)  ®C*)-1]  , 
where 

r " ( (vjkYj/k/)  j'  = 1 p ; k,k,=q2  + l q , 

and 

(4.3)  V (F)  = ( (o  . . , (F)  ) ) , 0 . . / (F)  = Cov(X  , . , X...) 

33  33  li^'i 
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Consider  now  a sequence  of  alternatives  fK  1 , where 

n' 


(4.4) 


K*  : 0 = (0,  n~\2) 


then,  • ~n(2)  * under  , has  the  same 

joint  distribution  as  that  of  S under  K*  . Noting  this 

— n n 

fact  and  using  the  results  of  Puri  and  Sen  (1969),  it  follows 
* 

that  under  has  asymptotically  a multinormal  distri- 

bution with  mean  vector  A[0,y2]C  = Ay[C(21) , C ] , and 
dispersion  matrix  A(F)®C  . Thus,  under  {1^}  , as  n -*  » , 

(4.5)  i(n_J5Sn(2))  -»  (AyC*  , A(F)®C*) 

Consequently 


(4.6) 

where 


r(vv 


pq- 


AC 


(4.7) 


&c  = Tr[f  • (T  (F)  ® C*)  "1] 


where  T (F)  is  given  by  (3.12). 

From  (4.2)  and  (4.7),  we  conclude  that  the  Pitman 

Asymptotic  Relative  Efficiency  (ARE)  of  £n  with  respect  to 

L is 
n 

(4.8)  L = AC/AL  = Trrf(T(F)  ®C*)_1lAr''r(v(F)  s'C*_1)l 
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_ * 

which  depends  on  r,F  and  C . If  F is  a mu Ltinorma 1 

cdf  and  if  we  use  the  normal  scores,  then  it  can  easily  bo 

checked  that  T ( F)  = 7(F)  and  hence  A = A . In  such  a 
~ CL 

case  the  normal  scores  test  and  the  normal  theory  likelihood 

ratio  tests  are  asymptotically  power  equivalent.  However, 

in  general  for  arbitrary  F , / is  bounded  by  the 

Z $ L 

minimum  and  maximum  characteristic  roots  of  V(F)T  1 (F)  , i.e. 

(4.9)  Ch  rv(F)T_1(F)-l  U t * Ch.  f V (F)  • T_1  (F)  1 

where  cIk  is  the  ithe  largest  characteristic  root.  (The 
bounds  of  v(F)T  ^(F)  may  be  studied  as  in  Sen  and  Puri 
(1967)  or  Puri  and  Sen  (1969).  Because  of  the  similarity  of 
the  work,  the  details  are  omitted) . In  passing  we  may 
also  remark  that  the  test  has  asymptotically  the  best 

average  power  with  respect  to  surfaces  in  the  parameter  space; 
it  has  also  asymptotically  the  best  constant  power  on  such 
surfaces  and  finally  it  is  asymptotically  most  stringed  test. 
The  proof  follows  as  in  Theorem  6.2  of  Furi  and  Sen  (l')(>9). 


5. 


ADF  Tests  for  parallelism  of 
(k) 

Let  , k = 1 , . . . , n^  be 


regression  sur faces . 
n^  independent  i vs 


with 


continuous  cdfs 
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(5.1)  Ffk)  (x)  = P[X.(k)  < x]  = F(x  - a,  - n c.(k)) 

1 ~ ~i  ~ --  -Tc  ~-k  -l 


1 < i < 


V 


k - 1 , . . . , s-  2 


We  desire  to  test  the  null  hypothesis 


(5.2)  H = 0 = ...  = g =g  (unknown) 

u — l — s — 


Here  the  0 's  are  pxt  matrices  and  the  c. 

— K — • i 


(k) 


are  t-vectors 


for  some  til  . A special  case  of^  p = t=l  has  been  studied 


in  detail  in  Sen  (1969).  If  we  let  g 


1 + _B  ^ * k = 1 , . . . , s , 


(so  that  0 ^ = 0)  , q = st,  then  the  result  follows  from  the 
theory  developed  in  section  3.  Therefore,  without  going  into 


the  details  of  derivation,  we  briefly  present  the  theory  here. 


For  the  kth  sample  [i.e.  X. 


(k) 


, i = 1, 


n^}  , define  the 


pxt  matrix  as  in  (2.2)  and  for  every  B f Rpt 


(k) 


S^'  (B)  as  in  ( 2 . 13 ) - (2 . 15)  . Let  then 


(5.3)  S (B)  = 

— n ~ 


5 

& 


(k) 


(B) 


n = / \ 

k£l  k 


Under  , we  estimate  the  common  B as  follows: 

(3.8)  and  (3.9),  we  let 


(5.4)  D = {b  • f V IS  . (b. 

~n  V j4i  8=1  n'jr 


) | = minj 


as  in 


(5.5)  g * center  of  gravity  of  D 
n J --  n 
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Let  then 


(5.6! 


; (k)  _(k)  . . . 

~n^  ~ Snk  (V  ' k = 1 s ' 


(5.7)  H(k)  = ( (S(k)  . S(k)  ,)  ) 

A V3*  V3  s “1 P ; r , r ' = 1 , 

(5.8)  c(k)=  y [c.(k)-  c l[c.(k)  -c  ]' 

“"k  fri  -1  - A ~x  A 


(5, 


(5  • 1Q)  C = Mn  ® Cj[k)  , k = 1 s . 

k ~k 

where  R(k)  is  the  rank  of  xfk)  - g . C_lk),...,g  . c(k) 

j.i  n,jlil  n,jt  it 

among  the  n^  aligned  observations  on  the  jth  variate  in  the 

kth  sample,  for  i-l,...,!^  ; j = l,...,p  ; k = 1 s . 

The  aligned  rank  order  test  statistic  for  testing  H in  (5.2) 


is  then 


(5.11) 


<■„  = ) TrrH(k)G-1l 


\ ~v 


Under  in  (5.2),  has  asymptotically  chi  square 

distribution  with  p(s-l)(t-l)  degrees  of  freedom  and  under 
the  sequence  of  alternatives  fK  } , where 


20. 


(5.12) 


c „(k) 


K : B V = # + n Y-u  . k = 1 s ; C 'I''  y =0 

n -k  x-  k , -nk  ~k 


k=l 


it  has  a non-centrality  chi  square  distribution  with 
p(s-l)(t-l)  degrees  of  freedom  and  non  centrality  parameter 


(5.13) 

where 


-1. 


= ) Tr[  I\  (T  (F)  ®C.  ) 
£ k^l  k ~ ~k 


(k)  (k) 


-1  (k) 


(5.14)  r = ( (y  . y 7 />).  l<k^S  and  C = Lim  ri  c 
k ]r  3 r ~k 

n-4cc 

which  we  assume  to  exist. 


\ 
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